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We reinvestigate the utility of time-independent constant mean curvature foliations for the nu- 
merical simulation of a single spherically-symmetric black hole. Each spacelike hypersurface of such 
a foliation is endowed with the same constant value of the trace of the extrinsic curvature tensor, K. 
Of the three families of isT-constant surfaces possible (classified according to their asymptotic be- 
haviors), we single out a sub-family of singularity-avoiding surfaces that may be particularly useful, 
and provide an analytic expression for the closest approach such surfaces make to the singularity. 
We then utilize a non-zero shift to yield families of A'-constant surfaces which (1) avoid the black 
hole singularity, and thus the need to excise the singularity, (2) are asymptotically null, aiding in 
gravity wave extraction, (3) cover the physically relevant part of the spacetime, (4) are well be- 
haved (regular) across the horizon, and (5) are static under evolution, and therefore have no "grid 
stretching/sucking" pathologies. Preliminary numerical runs demonstrate that we can stably evolve 
a single spherically-symmetric static black hole using this foliation. We wish to emphasize that this 
QO ' coordinatization produces .ft'-constant surfaces for a single black hole spacetime that are regular, 

static and stable throughout their evolution. 

> 
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I. CONSTANT CRUNCH SURFACES 

in 

■ In this paper, we address a single question: Is there a numerically-viable coordinatization of a Schwarzschild 
black hole spacetime foliated by hypersurfaces of constant (not necessarily zero) mean extrinsic curvature? In other 
words, can we coordinatize the Schwarzschild spacetime with constant mean extrinsic curvature (Tr(K) = constant) 
hypersurfaces so as to bound the growth of metric components and their gradients? We demonstrate here that the 
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I . single shift freedom yields a spacetime metric that is static, and therefore bounds the growth in time of such gradients. 
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A more complete analysis of the stability of our coordinatization, and a more thorough canvassing of the parameter 
space, will appear elsewhere, jjj Our foliation is consistent with that of Iriondo et al. [Q, who provided a generic 
constant mean curvature (CMC) foliation of the Reissner-Nordstrom spacetime for the purpose of finding trapped 
surfaces. In this paper we focus on the utility of CMC slicings for the numerical simulation of black holes, in support 
j_j ■ of the emerging field of gravity- wave astrophysics. 

The trace of the extrinsic curvature tensor (Tr(K) = K a a = K) at a point on a spacelike hypersurface measures 
the fractional rate of contraction of 3- volume along a unit normal to the surface. It represents the amount of "crunch" 
the 3-surface is experiencing at the point, at a given time. If all the observers throughout a spacelike hypersurface 
moving in time orthogonal to the surface experience the same amount of contraction per unit proper time, we say that 
the surface is a if-surface or a "constant crunch" surface. In this paper we examine foliations of a single spherically- 
symmetric, static black hole where each spacelike hypersurface has the same constant value of the extrinsic curvature, 
K. 

Generic if-surface foliations have found great utility in the numerical simulation of cosmological spacetimes. Q In 
addition to decoupling the three momentum constraint equations from the Hamiltonian constraint, these surfaces (in 
the case of compact or W- model universes) provide a convenient cosmological time parameter (K, or York, time). Q] 
Furthermore, for such cosmological spacetimes one has powerful existence and uniqueness theorems, p, u Extensive 
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work into the characteristics of these surfaces for Schwarzschild spacetimes has been done by Brill et al., M and 
foundational work into their use in numerical relativity was done by D. Eardley et al.. Q More recently Pervez et 
al. H provided a foliation partially covering the Schwarzschild spacetime with if-surfaces, with K ranging from — oo 
to oo, and Iriondo et al. provided a generic constant mean curvature foliation of the Reissner-Nordstrom spacetime 
for the purpose of finding trapped surfaces. In this paper we build upon the work of these investigators by examining 
the utility of these surfaces for numerical relativity in support of gravity- wave detectors. 

Although surfaces of constant K were thoroughly investigated decades ago, their use in current numerical simulations 
of black holes is conspicuously absent (apart from the use of maximal (K — 0) surfaces). jl(| One reason why 
these slicing methods have not been more fully developed is that they lag in time close in, to avoid crashing into 
the singularity, while they simultaneously evolve forward normally at the outer edge of the grid to allow for wave 
extraction. This tension, many fear, will unavoidingly lead to unbounded growth in the metric and extrinsic curvature 
components in the intermediate region, as is indeed found in maximal slicing. This computational concern has been 
referred to by the numerical relativity community as "grid stretching" or "grid sucking." We show in this paper that 
a proper choice of radial shift can yield a constant crunch foliation of a spherically-symmetric black hole without such 
pathologies. In fact, we foliate a Schwarzschild black hole such that the 3-metric and extrinsic curvature are both 
bounded and static (i.e. unchanging in time). 

To numerically evolve a black hole 3-space in time it is desirable to have a foliation, and its coordinatization, which 
satisfy the following four properties: 

1. Avoids black hole singularities or facilitates their excision. 

2. Possesses asymptotically null hypersurfaces to aid in radiation extraction. 

3. Minimizes steep gradients in the lapse, shift, 3-metric and extrinsic curvature tensor. 

4. Maximizes the future development of the initial data for the purpose of gravity-wave extraction. 

As a first step towards achieving these goals for systems containing multiple black holes, we explore the families of 
AT-surfaces in the Schwarzschild spacetime, and find a CMC foliation satisfying the above properties. 



In the next section we construct the fC-surfaces outside and inside the horizon. In section III we explore the 
properties of the if-surfaces, dwelling in particular on their approaches to the singularity. We also examine and 
illustrate the three families of if-surfaces. In section [fv| we derive a metric for Schwarzschild whose constant time 
slices are if -surfaces. We restrict our attention to a subfamily of if-surfaccs - surfaces which, when generalized to 
the colliding black hole spacetimes, support the gravity wave detection problem. We also present some preliminary 
numerical simulations using constant-crunch coordinates. We conclude with general comments on the applicability of 
if-surfaces to numerical calculations of more general black hole spacetimes. 



II. CONSTRUCTION OF CONSTANT CRUNCH SURFACES 

Over three decades ago Eardley and Smarr || carried out a generic classification of the spacetimes that could be 
simulated numerically, and investigated the limitations that the presence of singularities would impose. In their paper 
they argued that CMC slicings are particularly useful for numerical purposes. In particular, they demonstrated this 
explicitly by constructing numerical solutions to a wide array of dust collapse models. In a similar vein, Brill et al. § 
explored the nature of CMC slices of the Schwarzschild spacetime, and they also presented some numerical examples. 
In the present work, we explore the numerical utility of CMC slicings in the case of single black hole spacetimes. For 
the sake of clarity, we will commence with a re-derivation of the equations governing the CMC surfaces starting from 
Schwarzschild coordinates. From there, we will explore specific properties of the surfaces, paying particular attention 
to implications for numerical relativity. 

We wish to find a spacelike hypersurface in the Schwarzschild spacetime such that every point on the surface has 
the same constant value of the trace of the extrinsic curvature tensor. We have at our disposal the specification of 
the initial-value data, as well as the freedom to choose the lapse and shift throughout the evolution. To begin, let us 
take the standard coordinate system of a single black hole spacetime of mass M in Schwarzschild coordinates: 

ds 2 = -B{r)dt 2 + C(r)dr 2 + r 2 (d6 2 + sin 2 9 dip 2 ) , (1) 

with B(r) = (1 — 2M/r) and C(r) = l/B(r). It will be convenient to treat separately the regions inside and outside 
of the horizon. We will find that the two are related by an isometry. 
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r=R(t) 




r=2M 



n / t=T(r) 




FIG. 1: To construct CMC slices, at each point we introduce a unit normal vector n and the unit tangent vector £ to the 
surface. The local light cones are depicted by the light dashed lines. 



A. Outside the Horizon (r > 2M) 

Outside the horizon, CMC surfaces will be labeled by t — T(r) (Fig. |l|). The requirement that the trace of the 
extrinsic curvature be constant throughout this surface yields a first order differential equation for T(r), determined 
by examining the behavior of the normals to the surface. The normal n to the spacelike hypersurface T is given by, 



n = N V(t - T(r)) = n t dt + n r dr = N (dt - T'dr), 



(2) 



where No is a normalization constant and primes denote differentiation with respect to r. The normalization is fixed 
by demanding that 
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The contravariant components of the normal are given by 
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The trace of the extrinsic curvature is the fractional rate of contraction of 3-volume per unit proper time along the 
normal, namely 



K 
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r 2 dr 



(r 2 n r ) 



(11) 
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Substitution of n r into Eq. (11) yields a second-order ordinary differential equation for T. Integrating this equation, 
we find 



H = 



Br 2 T' \ 



y/C-BT 12 ) 

with H an integration constant and J an indefinite integral given by 



+ J, 



J = I KVBCr 2 dr = -ifr 3 . 



(12) 



(13) 



Along the surface the rate of change of proper time, dr, with proper distance, ds, is related to the slope of the surface 
T, 



dr B 



r. 



From Eq. (|T|) we find 



(dr 



{H-Jf 



\dsj (H-J) 2 + Br 



(14) 



(15) 



B. Inside the Horizon (r < 2M) 

Finding the ^-constant slices of Eq. (Q) within the horizon is similar to the calculation done in the previous section; 
however, as the roles of time and space coordinates reverse within the horizon, we will find it useful to parameterize 
our spacelike surface as a function of coordinate t, and look for if-constant surfaces of the form (Fig. |l|) 

r = R(t). (16) 

The normal n to the spacelike hypersurface R is given by 

n = N V{R(t) - r) = n t dt + n r dr = N (Rdt - dr), (17) 

with differentiation with respect to t denoted by dots and the Nq a normalization constant fixed by: 

n-n = -1 (18) 
= g u n t n t + g rr n r n r (19) 

->*(h-Ti#)- (20) 



We have, therefore, 



and 



^Vo = , . (21) 
VCR 2 - B 



(22) 



Vcr 2 - b 



Vcr 2 - b 

The contravariant components of the normal are given by 



n t = ; = ■ (23) 



n r = g rr n r = 1 (24) 

CVCR 2 - B 

n l = g u n t = R (25) 

BVCR 2 - B 
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From Eq. (jllj), we once again find that fixing the trace of the extrinsic curvature gives us a second-order differential 
equation for R(t), namely, 

o 2CB'R 2 -b(b' + C'R 2 + 2CR] 

crVcrz-b 2(ciV-B)' 

which can be simplified to 

KR*R=-±( - ). (27) 

dt \Vcr 2 -bJ 

Paralleling the approach from the last section, we introduce an integration constant, H, and an indefinite integral 
J (given by Eq. (13)), to obtain the first integral: 



H = = + J. (28) 

V Vcr 2 -b 



From Eq. (g7|) we find that the "proper velocity" along the surface, ds/dr = y/C/BR, results in the same equation 
both inside and outside of the horizon (Eq. (pi])). This can be rewritten as 

ds\ 2 BR 4 

^)= 1+ (tf-7F (29) 

The spacelike if-surfaces obtained from the first integrals, Eqs. ( |l5| ) and (29), differ only by an isometry, 

T' I. (30) 



III. PROPERTIES OF THE A-SURFACES 



The spatial metric of a _ftT-surface outside of the horizon is given by 



ds 2 = d£ 2 +r 2 dn 2 (31) 

= (C - BT' 2 )dr 2 +r 2 dn 2 . (32) 

Within the horizon it becomes 

ds 2 = d£ 2 + r 2 dn 2 (33) 

= (CR 2 - B\ dt 2 + r 2 dn 2 . (34) 



These two expressions differ by the isometry of Eq. (|3C|). Using Eq. (15) we can rewrite them in terms of H and K 



r 
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ds* = n dr z + r A d£l\ (35) 

(H - J) 2 + Br* 

From this we arrive at the scalar curvature of the if-surface: 

<3>fl + (36) 
3 r° 

Similarly, by using Eqs. ( |l5| ) and (^9j), the extrinsic curvature associated with observers moving on world lines 
orthogonal to the iT-slices are also expressible in terms of K and H: 

Kt = \ K + 2 J, (37) 

4 = K\ = U-^. (38) 
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Schwarzschild Kruskal-Szekeres 
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FIG. 2: An example of the three families of spacelike if-surfaces for K — —0.1. The first row (HH) depicts a representative 
horizon-to-horizon surface using H = —1.25, which corresponds to R m in ~ 1.816 . This spacelike hypersurface is represented 
both in Schwarzschild coordinates (left column) and in Kruskal-Szekeres coordinates (right column). The middle two graphs 
are the horizon-to-singularity (HS) surfaces using H — —1.43. The bottom two graphs represents a typical singularity-to- 
singularity (SS) surface. We have used H = —1.25 to generate this SS /('-surface. The HH and SS surfaces are close to their 
critical radii (R c ~ 1.5646) and therefore appear flattened, as described in the text. 



The -ftf-surfaces are therefore parametrized by two constants: the trace of the extrinsic curvature tensor, K, and the 
constant of integration, H. In addition one must fix a single point on the surface, t Q — T(r a ), which amounts to 
setting a time translation parameter. As can be seen in Eqs. (p6|)-(^8|), the constant H controls the variation of the 
intrinsic and extrinsic curvatures over the ^-surface. 

To elucidate the nature of the ^-surfaces, we numerically integrate Eqs. (|l5|) and (|29|). We find that within the 
horizon there are 3 classes of if-constant surfaces, differentiated by their asymptotic behavior. The singularity- 
singularity surfaces (SS) begin at the singularity aligned with the null surface, reach up towards the horizon, and 
then fall back, reaching the singularity along the null cone. The horizon-horizon surfaces (HH), which we have also 
dubbed "horizon- hugging" surfaces, asymptote to the horizon (r — > 2 M for |t| — >■ cxd in Schwarzschild), dipping down 
towards the singularity in between. This feature was previously remarked upon by Brill et al. M. The asymptotes 
converge toward a null surface at the horizon. Finally, the horizon-singularity (HS) surfaces begin at the horizon, and 
asymptote in to the singularity. Representative surfaces for the value K = —0.1 are shown in Fig. |[ We integrate 
the HH and HS surfaces across the horizon into the region r > 2M by imposing continuity of the surface and its first 
derivative at the horizon. Because of the isometry, Eq. (|30|), the surfaces outside of the horizon are characteristically 
similar to those on the inside; in particular, both sets are null at their asymptotes. 

We have chosen to use the acronym HH for the horizon-to-horizon hypersurfaces, in lieu of referring to them as 
"regular" hypersurfaces Q|, as each of the three types of fC-surfaces are, in a strict sense, regular. In particular, each 
surface asymptotes to a null surface, be it at the singularity or the horizon. Observers on such a surface, or more 
precisely, observers that are time synchronized throughout the surface, are never seen crossing the horizon, nor do 
they ever reach the singularity! Outside the horizon, every HH and HS K-surface (K ^ 0) asymptotes for large r 
to null infinity. X-surfaces corresponding to positive values of K asymptote to past null infinity, and asymptote to 
future null infinity for K < 0. 

To gain a qualitative understanding of the X-constant foliation, it is useful to analyze Eq. (|2^ ) as an energy 
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FIG. 3: Phase space diagram of the classification of K -surfaces. The Tf-surfaces have been classified into 3 groups according 
to their asymptotic behavior (see text). Which class a particular surface falls into depends upon the values of K and H, as 
detailed in the figure. For the purposes of numerics, the surfaces of interest are to be found in the "Future HH and SS" wedge, 
with K < (see section |tl| ). 



conservation equation for a particle of unit total energy (E = 1) moving in the potential 

(1 



V(r) 



2M\ r 4 



{H-\Kr*) 



2 ■ 



(39) 



By using this energy equation we can determine the closest approach to the singularity, R m in, of a given HH in- 
constant surface. Two conditions must be satisfied to determine R m in- First, the closest approach occurs when 
R = (ds/dr) = 0, which is equivalent to demanding 



V(R m in) — 1 

This condition leads to a sixth-order polynomial in R m m- 

K 2 



- R 4- R 

min * min 



H = 0. 



(40) 



(41) 



Second, the solution for this surface at R(t m i n ) = Rmin must be concave, so as to rule out the SS surfaces (which 
bend towards the singularity rather than the horizon). This is enforced by demanding: 



R(U, 



(R n 



2) 



-3 + 2R mm + KRl 



\R(t„ 



R 3 



> 0. 



(42) 



The solution of these two conditions, as shown in Fig. ^, gives rise to the emergence of two critical values for H, 
namely H±, for a given value of K. In addition, an HH surface can be made to approach arbitrarily closely to the 
singularity at r = by choosing an appropriately large positive value of K. Negative values of K tend to "hug the 
horizon." 

For fixed K and H we know how to compute how closely a CMC surface comes to the singularity. But, for a given 
value of K, what value of H gives the closest overall approach to the singularity? We can determine the critical values 
for R and H, given by R c and H± respectively, by looking at the point where the first and second time derivatives 
of R(t) vanish. ii_ occurs along the lower boundary of the contour plot in Fig. and corresponds to the _ftT-surface 
that reaches down the furthest towards the singularity for a given value of K . The vanishing of R\ p_ leads to the 
following equation for R c : 



(R c - 2) -3 + 2R 



K! " :aI ~r~ 1 1 = °" 



(43) 
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FIG. 4: Critical value of R c = R- as a function of K for Schwarzschild coordinates and M = 1. One can readily see that the 
Tf-surface "hugs" the horizon for large negative values of K. 



which can be rewritten as a 4th order polynomial in R c 

K 2 R 4 r - 2K 2 Rl + 4i? 2 - 12R C + 9 = 0. 



(44) 



One must take care in examining the roots of this equation, as there are more solutions to Eq. ( fj4| ) than there are for 
Eq. (f43|). Nevertheless this equation gives two distinct real roots, depending on the sign of K: 



R 



c\K>0 



16 



__ Vx 



i 1 1 / 8 2 16 y/x 

R c \k<o - - + -J3-^- X + — + -^= + — , 



where 



£ ee 32 + 108Jf 2 + 243X 4 + 21K 2 y/lQ + 56K 2 + 81X 4 , 
8 16 + 36iv 2 2 1 / 3 £ 1 / 3 



X 
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3K 2 3 2 1 /3iT 2 ,ci/3 ^ 3K 2 



(45) 
(46) 

(47) 
(48) 



When K = we see from Eq. ( fl3|) that R c = 3/2. The regular if -surfaces are thus bounded between the horizon at 
R + = 2M and R- = R C M. Using Eq. (p8|), and setting 7? = 0, we obtain, for the case of a black hole in Schwarzschild 
coordinates 



H± = 



B±R 2 ± KR 3 ± 



+ 



with 



We therefore have 



2M 

= ~R± 



-M 3 K. 



(49) 



(50) 



(51) 



For large values of \K\ one can show that the critical value of R, depends upon the sign of K . In particular, 

n — : ;„.r.i > , .. ... . (52) 



2 - |X- 2 for K < -1 
(|) 1/3 K~ 2 / 3 forK»l 



This in turn gives the following asymptotic values for 

-I 
2if 



2K 



for Jf < -1 
for K > 1 



(53) 
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IV. CONSTANT CRUNCH COORDINATES: A SPACETIME METRIC FOR A If-SURFACE 
FOLIATION OF THE SCHWARZSCHILD BLACK HOLE. 



A number of features of the if-constant surfaces presented in the previous sections seem particularly well suited 
to the numerical analysis of generic black hole spacetimes. First, the surfaces asymptote to a null surface, making 
them effective for gravity wave extraction. Second, the if -surfaces naturally avoid the crushing singularity. Finally, 
for large negative values of if the surfaces "hug the horizon." This last feature, illustrated in Fig. [|, allows one to 
focus attention on the region relevant for gravity wave generation — the region outside the horizon. 

In this section we generate a if -constant foliation for the Schwarzschild black hole that, in addition to the properties 
just mentioned, also has regular and static metric and extrinsic curvature components. To generate this if -constant 
slicing we use the coordinate transformation 

i = t-T(r), (54) 
p = r. (55) 

Under this transformation, the metric from Eq. (hi) becomes 



ds 2 = _ (1 _ ^_ )d p o (J H) p4 

P ! / (ff _J)2 + (1 _2 M) ^4 (tf-J)2 + (l_2M )p 4 / • 



P 



The constant t slices of this metric are if -constant surfaces. It is to be noted that Eq. (|56|) agrees with Eq. (53) of 
Iriondo et al. (for the case of constant if and vacuum). However, in order to regularize the g pp metric component 



at the throat, we add the isotropic- like radial transformation |11 

1/1 



, M ^Rmin + P + \Z-RrainP + P 2 j , (57) 

with Rmin the minimum coordinate location of the throat, given by Eq. (fll]). This coordinate representation of a 
black hole spacetime provides a foliation with if-constant, ii-constant spacelike hypersurfaces. Each hypersurface is 
metrically equivalent to all others — the surfaces are independent of t, and hence static. In addition, the hypersurfaces 
are asymptotically null (T'(r) — + 1 as r — * oo). Furthermore, the lapse, the shift, and all of the 3-metric and extrinsic 
curvature components are regular and well behaved, as illustrated in Fig. pi 

In addition to restricting ourselves to the singularity-avoiding family (HH) of if -surfaces, two additional conditions 
on the if -surfaces are demanded by the nature of our problem - the eventual simulation of the gravity-wave emission 
from two interacting black holes. First, the foliations must asymptote at large f to future null infinity. Therefore, we 
must restrict our attention to negative values for the trace of the extrinsic curvature tensor, if. Second, we require 
that such negative if hypersurfaces enter the future singularity region. This will ensure proper coverage of the relevant 
region just above the future horizon, which is precisely where the gravity waves are produced. However, we expect the 
initial-data formulation for such surfaces to be involved, and this may guide our choices even more systemically. The 
two additional requirements limit us to the relatively narrow wedge of Fig. formed by restricting to the "Future 
HH and 55" shaded region with if < 0. 

A representative constant-crunch foliation generated by Eq. (|56| ) is shown in Fig. ^[ The avoidance of "grid 
stretching" is accomplished by a suitable choice of shift vector. To illustrate the non-zero shift we show the if = — 1, 
H w —3.11 foliation of Schwarzschild in Fig. with the explicit misalignment of the r = constant line segment and 
the normal vector. 

Finally, we describe several preliminary numerical experiments using if -surfaces. A full treatment of a single 
black hole using if -constant foliations will be presented elsewhere Here we present several sample evolutions 
demonstrating the utility of constant mean curvature slicings. Figures || and ^ display results from the simplest 
possible test of if -constant foliation of the Schwarzschild geometry; the domain is taken to be a thin shell close to 
the horizon (in this case r € [1,5]), analytic Dirichlet conditions are applied at both boundaries of the computational 
domain, and analytic lapse and shift conditions obtained from Eqs. (56) and (|5^) are used. The figures represent 
the singularity- avoiding foliation K = — 1 and H = — 3, for which the evolution was found to be stable and accurate 
over very long timescales. Using 50 grid points, the code succesfully ran beyond t — 50, 000M while maintaining high 
accuracy. The fractional error in the metric components was typically 1 — 2%. 

Fig. g shows the convergence of the mean fractional error in the metric component a — ^/<?ff as a function of time. 
Each curve has been rescaled by a factor of 4 P , where the number of grid points is given by 400/2 p , with p = 0, 1, 2, 3. 
For t > 2 the solution is approximately second order accurate. Initially, noise generated on the inner boundary causes 
fluctuations whose magnitude is largely independent of the number of grid points. Fig. ^| shows snapshots of the 
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FIG. 5: The radial behavior of the various components of the spacetime metric (Eq. (pqjj, under the isotropic-like transformation 
given in Eq. for M = 1, K = —0.1, H = —1.0 and _R m i n ~ 0.47922. From upper left to lower right we plot the radial 

metric component (gff), lapse (a), radial shift (/3 r ), diagonal components of the extrinsic curvature tensor {K^ and Kg = K^), 

and in the lower right frame the trace of the extrinsic curvature tensor (K = K? + Kg + K^) as a consistency check. All of 
the functions are regular and well behaved. The growth of the lapse and shift for large f is expected, as the surfaces become 
asymptotically null. We show an exploded view of the behavior of the radial shift near the throat to emphasize that the shift 
changes sign and becomes positive before reaching f — R m i n - 



Schwarzschild Kruskal-Szekeres 
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FIG. 6: The constant crunch coordinate foliation of the Schwarzschild spacetime. We show the foliation generated by Eq. (|H(j) 
for M — 1, K — —0.1, H = —1.25 and Rmin ~ 1.816, in Schwarzschild and Kruskal-Szekeres coordinates. 

Hamiltonian constraint at various times in the evolution. Qualitatively, a wave which is triggered by truncation error 
is seen to propagate outwards from the inner boundary. The amplitude reduces rapidly, before growing once more 
as it is reflected off the outer boundary. The magnitude and speed of propagation of the wave quickly decay as the 
wave moves back into the domain, leaving a static solution which is stable beyond t = 50, 000M. The numerical error 
which initially propagates through the domain is caused entirely by the analytic Dirichlet boundary conditions, and 
can be largely eliminated by the use of more realistic conditions. |jj 

The numerical runs presented here evolve a K = — 1 singularity-avoiding hypcrsurface that asymptotes to null 
infinity, entering the past singularity region. This is not of the class of if -surfaces emphasized in this paper for use in 
numerical relativity. Ideally, we would have preferred presenting the evolution of a K = —1 HH surface that enters 
the future singularity region. However, we were unable to find an HH surface in the future singularity region that 
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Schwarzs child 




FIG. 7: The non-zero shift in the constant- crunch coordinatization of the Schwarzschild spacetime. We show a foliation generated 
by Eq. (Q) for M = 1, K = — 1, ~ —3.11, and R m in ~ 1.9. We show the r = constant line segment at a point on one of the 
surfaces, together with the inward-pointing normal vector. The misalignment of these two line segments indicates the non-zero 
shift in this static coordinatization. A magnified view of the two line segments is shown in the bottom panel. 




FIG. 8: A numerical example: convergence of the radial metric component for the evolution of an HH surface spanning the 
past singularity region and reaching positive null infinity. We show the time evolution of a single black hole with M = 1, K = 
— 1,H = —3, and plot the mean fractional error in the radial metric component a = ^/gTf- Results are shown for four different 
resolutions, iV = 400/2 p (p = 0, 1, 2, 3), where each curve has been rescaled by a factor of 4 P . For t > 2 the solution displays 
approximately second order convergence, with a fractional error of around 0.1% for N = 200. Long term evolution is stable, 
and the system has been succesfully evolved beyond t = 50, 000M. 
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FIG. 9: Snapshots of the Hamiltonian constraint. For the same parameter choices as in Fig. ^ we display the Hamiltonian 
constraint during the early stages of evolution. A wave generated by truncation error at the inner boundary propogates 
outwards, is reflected by the outer boundary back into the domain, and then rapidly decays. No further evolution is seen in 
the Hamiltonian constraint after t « 5. 



evolved stably using the naive Dirichlet boundary conditions and analytic lapse and shift conditions. Nevertheless, we 
have evolved such surfaces stably by incorporating (1) area locking shift conditions, (2) isometery conditions at the 
throat, and (3) an outgoing boundary condition based on the difference between the numeric and analytic solutions. 
These numerical results, and the corresponding stability analysis, are not within the scope of this paper and will be 
presented elsewhere, pf 



V. FROM ONE BLACK HOLE TO TWO SPINNING BLACK HOLES 



In this paper we considered a static spacetime metric for the Schwarzschild black hole, with spacelike hypersurfaces 
of constant (not necessarily zero) value of the trace of the extrinsic curvature tensor. This slicing provides a natural 
generalization of the maximal slicing scheme currently in use in some numerical approaches to the binary black hole 
problem. 

An essential feature of our if-constant metric is a spatially-varying radial shift vector, which allows the surfaces 
to avoid the singularity while evading the grid stretching problems often encountered with other metrics. Work is in 
progress to develop a geometric handle on our shift condition, akin to that of the "minimal distortion" shift often 
discussed. The inner and outer boundary conditions are also particularly convenient. In the inner regions, the "horizon 
hugging" feature of the if-slices, together with their regularity, may remove the need to excise the grid within the 
apparent horizon, thus providing a natural "boundary without a boundary" avoidance of the singularity. In addition, 
at the outer boundary the surfaces are asymptotically null, which may aid in gravity wave extraction. 

In examining the characteristics of the if-slices we reviewed the three families of -fT-surfaces, including the horizon- 
to-singularity (HS) surfaces, as well as the more familiar horizon-to-horizon {HH) and singularity-to-singularity (SS) 
surfaces. All three families of surfaces either asymptote to the singularity, or to spatial infinity along a null surface. 
Thus the family of observers that make up such a foliation never "observe" the surfaces reach the singularity. This 



13 



suggests a possible payoff in using the HS surfaces in the evolution of black hole spacetimes, as there is no danger of 
the surface reaching the singularity. Just as there are natural boundary conditions which "freeze out" gravitational 
radiation as it progresses to null infinity, so too can we expect boundary conditions where the dynamic space freezes 
out as it approaches the singularity. Are there such natural boundary conditions at the future singularity? It would 
be worthwhile to investigate such asymptotically-null boundary conditions at the singularity of an HS surface. 

The future utility of numerical simulations of black hole spacetimes hinges in large part on a suitable choice of 
coordinates for the initial data, and on the particular evolution of these coordinates through the four lapse and 
shift conditions. These conditions are the only handles by which to manage the growth of the metric and curvature 
components during evolution. It is through judicious choices of lapse and shift that one is able to effectively enable 
singularity avoidance, and allow for efficient extraction of gravitational radiation. The simplistic example presented 
here may provide some guidance as to how to proceed in the more general case. A"-surfaces in the generic two black 
hole problem can be expected to preserve the singularity-avoidance "horizon hugging" behavior, as well as remain 
asymptotically null at the boundaries. It remains to be seen what shift vectors are required to manage the growth of 
the intrinsic and extrinsic properties of the metric in the more general case, though the constant-crunch shift presented 
here is a natural starting point. 

We are currently pursuing three avenues towards furthering and generalizing the work presented here. First, we are 
analyzing the stability of this coordinatization to small perturbations. Second, we are using the (l-fl)-dimensional 
code discussed in the last section to examine a greater portion of the parameter space of initial data, so as to fully 
explore the numerical stability of the slicing. From our preliminary numerical investigations, we expect to be able 
to do a full spacetime "evolution" of Schwarzschild, and have it run stably and accurately for extended periods of 
time over a wide range of H-K parameter space. Finally, we are analyzing the Oppenheimer-Snyder collapse in such 
A-slicings, following the lead of Eardley and Smarr, || as this will further test our slicing in a non static setting. Q 
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